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摘 AAA 精确 度 不 高 、 收 敛 速度 慢 等 问题 进行 了 研究 ， 提 
出 一 种 多 模 态 多 目标 差分 进化 算法 。 首 先 将 非 线性 方程 组 转换 为 多 模 态 多 目标 优化 问题 ， 初 始 化 一 个 随机 种 群 并 对 种 
ey 然后 通过 非 支配 解 排序 和 决策 空间 拥挤 距离 选择 机 制 ,挑选 种 群 中 的 一 半 优 质 个 体 进 行 变 异 ; 
接着 在 变异 过 程 中 采用 一 种 新 的 变异 策略 和 边界 处 理 方 法 以 增加 解 的 多 样 性 ; 最 后 通过 交叉 和 选择 机 制 使 优质 个 体 进 
行进 化 ， 直 到 搜索 到 全 部 最 优 解 。 在 所 选 测 试 函 数 集 和 工程 实例 上 的 实验 结果 表明 ， 该 算法 能 有 效 地 搜索 到 非 线 性 方 
程 组 的 解 ， 并 通过 与 当前 四 个 算法 进行 比较 ， 该 算法 在 解 的 数量 和 成 功率 上 具有 优越 性 
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Multimodal multi-objective differential evolution algorithm for solving nonlinear equations 
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Engineering, Zhongyuan University of Technology, Zhengzhou 450007, China) 


Abstract: Aiming at the problems such as incomplete solutions, low accuracy and slow convergence speed of current algorithms 
in solving nonlinear equations, this paper proposed a multimodal multi-objective differential evolution algorithm. Firstly, it 
transformed the nonlinear equations into multimodal multi-objective optimization problems, and initialized a random population 
and evaluated all individuals in the population. Through the non-dominated sorting and decision space crowding distance 
selection mechanism, it selected half of the individuals in the population to mutate. Then it used a new mutation strategy and 
boundary processing method to increase the diversity of solutions. Finally, it evolved the high-quality individuals through 
crossover and selection mechanism until it found all the optimal solutions. The experimental results on selected test function 
sets and engineering examples show that the algorithm can effectively search for optimal solutions and it is superior to the other 
four algorithms in the number of solutions and success rate. 
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生 。 与 传统 数学 方法 相 比 ， 进 化 算法 是 一 种 成 熟 的 


鲁 棒 性 和 广泛 适用 性 的 全 局 优化 方法 ， 具 有 自 组 织 、 自 适应 、 


Fd 


直 以 来 都 是 数学 学 科 中 的 一 个 难题 。 ”自学 习 的 特性 ， 能 够 不 受 问题 性 质 的 限于 


， 有 效 地 处 理 传统 数 


在 实际 生活 中 ， 很 多 物理 、 生 物 、 工 业 里 复杂 的 问题 都 可 以 看 。 学 方法 难以 解决 的 复杂 问题 。 因 此 ， 利 用 进化 算法 求解 非 线 性 
改 是 非 线性 方程 组 问题 ”“”。 考 虑 到 非 线性 方程 组 本 身 的 复杂 性 ，，” 方程 组 已 成 为 当前 一 个 重要 的 研究 方向 。 目 前 应 用 于 求解 非 线 
如 何 有 效 地 获得 非 线性 方程 组 的 所 有 解 已 成 为 目前 研究 的 主要 ”性 方程 组 的 主要 进化 算法 有 遗传 算法 ”、 差 分 进化 算法 "“、 人 工 
目标 。 传 统 数学 方法 如 牛顿 法 、 拟 牛顿 法 、 欧 拉 法 、 复 合 形 法 、 鱼 群 算法 ”、 人 工蜂 群 算法 ”、 人 工 董 火 虫 算 法 ”“、 布 谷 鸟 算法 
梯度 法 、 和 迭代 法 、 共 斩 方 向 法 、 混 沌 法 、 最 速 下降 法 、 单 纯 型 ” “等 。 

法 等 ， 大 多 对 问题 的 初始 值 具有 很 严格 的 要 求 ， 另外 对 于 某 些 差分 进化 算法 作为 当前 流行 的 进化 算法 之 一 ， 由 于 其 容易 
导数 不 存在 或 是 导数 难 求 的 方程 ， 传 统 数 学 方法 也 具有 一 定 的 。 实现、 可 控 参数 少 等 特点 ， 已 经 被 一 些 研究 学 者 用 来 求解 非 线 
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性 方程 组 。 陈 海 霞 等 人 ”提出 将 非 线性 方程 组 转换 为 一 个 不 可 
微 优 化 问题 ， 并 用 光滑 函数 代替 不 可 微 的 极 大 值 函数 ， 最 终结 
果 表 明 该 算法 在 求解 非 线性 方程 组 的 准确 性 和 有 效 性 上 均 优 于 
其 他 同期 算法 。 针 对 差分 进化 算法 进化 后 期 收敛 缓慢 和 稳定 性 
不 强 的 缺陷 ， 刘 利 斌 等 人 "将 拟 牛 顿 法 插入 到 差分 进化 算法 当 


该 


AN 
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iv | 
许 伟 伟 ， 等 :多 模 态 多 目标 吝 分 进化 算法 求解 非 线 性 


min g(x) A We A | 


min 8 (0)= | (0) 


min g,(x) i | 


方法 将 方程 组 中 每 一 个 方程 都 当成 一 个 优化 目标 ， 


中 ， 有 效 提高 了 差分 进化 算法 的 收敛 速度 和 重 棒 性 。 封 全 喜 等 ”方便 。 
人 ”通过 代理 约束 概念 和 修正 极 大 粒 函 数 ， 将 非 线 性 方程 组 转 该 方法 
换 为 无 约束 优化 函数 并 利用 自 适 应 正 交 差分 进化 算法 求解 ， 验 。 ”方法 是 


证 了 所 提 算 法 的 有 效 性 。 
上 述 算 法 在 一 定 程度 上 实现 了 非 线 性 方程 组 的 求解 ， 但 仍 
然 存 在 所 得 解 的 个 数 不 完 整 、 精 确 度 不 高 收敛 速度 慢 等 问题 。 
为 解决 这 些 问题 ,本 文 提出 了 一 种 多 模 态 多 目标 差分 进化 算法 。 
首先 将 种 群 初始 化 并 计算 每 个 个 体 的 适应 度 值 ， 然 后 通过 一 种 
预选 择机 制 将 优质 个 体 筛选 进行 变异 ， 接 着 在 变异 过 程 中 采用 
新 的 变异 策略 和 边界 处 理 办 法 提高 搜索 解 的 多 样 性 ， 最 后 经 过 


但 随 着 方程 个 数 的 增多 ， 求 解 的 复杂 度 将 会 增 大 。 


该 
组 的 解 
择机 制 
标 转 


交叉 选择 优质 个 体 进行 多 次 进化 ， 最 终 搜索 到 全 部 最 优 解 。 实 
验 结果 表明 ， 该 算法 能 够 有 效 提 高 非 线 性 方程 组 的 解 的 数量 和 


1 问题 描述 和 转换 方法 


一 般 地 ， 一 个 非 线性 方程 组 可 以 被 描述 如 下 : 


BR ee eh a! 


fsXa, Nase Nh) =0 


(1) 


fn Ks Ns hh) =0 


其 中 : m 代表 方程 的 个 数 ，n 代表 决策 变量 的 个 数 。 
进化 算法 在 求解 非 线性 方程 组 时 需要 将 非 线 性 方程 组 转换 


为 优化 问题 。 当 前 转换 的 方法 主要 单 目标 优化 、 约 束 优化 和 多 
目标 优化 有 三 种 。 第 一 种 是 将 非 线性 方程 组 转换 为 单 目标 优化 
问题 ， 它 的 形式 如 下 : 


min P(x)= MG (2) 
声 
或 者 
min P(X)= Ff20, x, x. 7x0) G) 
A 


第 二 种 是 将 非 线性 方程 组 转换 为 约束 优化 问题 ， 它 的 形式 
如 下 : 


min PO)= Df )| (4) 


subject to fi(X,N, Xa Kh,) 0,17=1,2,...,m 


这 两 种 转换 方法 的 共同 点 是 都 有 一 个 优化 目标 ， 得 出 的 结 


该 


一 个 权 值 ， 能 


值 双 


双 目 标 转换 法 ， 它 的 表达 式 如 下 : 


min 万 (2) = 为 + 六 

min f(x)=1—an | 
方法 将 线性 方程 与 单 目标 优化 方程 结合 ， 
转换 成 线性 Pareto 前 沿 上 的 解 ， 直 观 
的 作 上 
换 法 ， 


它 的 表达 式 如 下 : 


>or ， 
局 | 
i=l 


0 


min fi(x)= 


i=1 
OX 
min (0)=1- Eo— + Df ,No,) 
@, i=1 
i=1 


方法 在 双 


(5) 


简单 
丸 此 


只 适用 于 具有 少量 方程 的 非 线性 方程 组 的 求解 。 第 二 类 


(6) 


把 非 线性 方程 
9 效 。 但 是 由 于 选 
容易 造成 部 分 解 的 丢失 。 第 三 类 方法 是 多 权 值 双 


(7) 


标 转换 法 的 基础 上 对 每 一 个 未 知 变量 都 设 定 


了 效 避 免 最 优 解 的 丢失 。 因 此 ， 本 文采 


多 权 


标 转换 法 求解 非 线 性 方程 组 。 


2 ”标准 差分 进化 算法 


差分 进化 算法 是 一 门 新 兴 的 进化 计算 技术 ， 主 要 用 于 求解 
连续 变量 的 全 局 优化 问题 。 该 算法 通过 变异 、 交 叉 、 选 择 三 种 
操作 不 断 地 和 迭代 计算 ， 保 留 优良 个 体 ， 淘 汰 劣质 个 体 ， 直 到 找 
到 全 局 最 优 解 。 差 分 进化 算法 具有 操作 简单 、 控 制 参数 少 、 易 


编程 实 
以 及 非 


现 等 优点 ， 同 时 它 又 尤其 适合 求解 多 变量 、 非 凸 、 
线性 优化 问题 。 其 基本 的 算法 流程 如 下 : 


多 峰 


a) 初 始 化 种 群 个 体 和 相关 参数 。 主 要 包括 种 群 大 小 、 交 叉 


率 、 缩 放 因 子 、 最 大 评价 次 数 和 终止 条 件 等 。 
bb) 变异。 通过 变异 策略 将 当前 个 体 进 行 变异 。 有 具体 操作 是 


通过 把 种 群 中 两 个 成 员 之 间 的 加 权 差 向 量 加 到 第 三 个 成 员 上 ， 


以 此 来 产生 新 参数 向 量 ， 从 而 产生 新 个 体 。 


9 交叉 。 变 异 向 量 的 参数 与 另外 预先 确定 的 目标 向 量 


按 一 定 


参数 


规则 混合 来 产生 实验 向 量 。 
dd) 选择 。 若 实验 向 量 的 适应 度 值 比 目 标 向 量 的 适应 度 值 好 ， 
量 束 在 下 一 代 中 代替 目标 向 量 ， 否则 ， 保 留 目 标 向 量 。 


实验 向 


果 是 一 个 最 优 解 。 而 要 想 同时 找到 多 个 最 优 解 ， 需 要 采用 第 三 
种 转换 方法 : 多 目标 优化 。 当 前 将 非 线 性 方程 组 转换 为 多 目标 
优化 问题 的 形式 主要 有 三 类 ， 第 一 类 是 直接 转换 法 ， 它 的 表达 
式 如 下 : 


对 每 


个 体 都 进行 选择 以 此 得 到 新 种 群 。 


6) 循环 。 直 到 满足 终止 条 件 ， 停 止 。 


3 ”多 模 态 多 目标 差分 进化 算法 


多 


模 态 多 


标 优 换 问题 ， 是 指 在 决策 空间 中 至 少 有 两 个 
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Pareto 解 集 对 应 


目标 空 


间 上 的 同一 个 Pareto 前 沿 . 如 图 1 所 示 ， 


决策 空间 上 的 Pareto 解 集 1 和 Pareto 解 集 2 同时 对 应 目标 空间 


上 的 同一 条 Pareto 前 沿 。 
22| ”Pareto 解 集 1 


广 


Pareto 前 沿 


Chi 
ee 


许 伟 伟 ， 等 
的 整数 ， 且 记 天 户头 太 关 凡凡 i。 纵 放 医 
非 支配 解 排序 拥挤 距离 选择 


(目标 空间 ) 


Ey 


Fl 
F2 


(决策 空间 ) 


-Et 
Wg 


合作 其 


解 非 线性 


乔 


方程 组 


子 F 的 取 值 是 [0,1] 。 


Pareto 解 集 2 P P 
决策 空间 标 空间 图 2 个 体 预 选择 机 制 
图 1 多 模 态 多 目标 问题 . . 
就 传统 差分 进化 算法 而 言 ， 对 于 超出 了 搜索 边界 的 个 体 往 
而 非 线性 方程 组 中 也 有 一 些 解 的 分 布 情况 与 多 模 态 多 目标 往 将 其 置 为 边界 值 。 而 对 于 MMODE， 该 算法 加 入 了 一 种 新 的 
问题 中 的 情况 类 似 。 基 于 这 种 思想 ， 本 文 提出 了 一 种 多 模 态 多 ”边界 处 理 方法 。 它 将 停留 在 边界 上 的 个 体重 新 分 布 在 整个 决策 
目标 差分 进化 算法 (MMODE)。MMODE 在 标准 差分 进化 算法 ”空间 中 ， 目 的 是 在 搜索 过 程 中 保证 解 的 多 样 性 。 这 种 新 的 变异 
基础 上 又 增加 了 一 种 个 体 预选 择机 制 和 一 个 变异 边界 处 理 方法 ， ”策略 表达 式 如 下 : 
这 样 可 以 有 效 提 高 搜索 精度 和 效率 。 它 的 具体 步骤 如 下 ; Vo = Xue -PO Xs) -F*CK Na) (9) 
a) 初 始 化 。 生 成 一 个 种 群 大 小 为 w 的 随机 种 群 P， 将 缩放 车 采用 这 种 变异 策略 变异 后 新 个 体 仍 然 超出 边界 ， 就 将 该 
因子 、 交 叉 率 、 最 大 评价 次 数 等 赋 初 始 值 。 个 体 随 机 分 布 在 种 群 中 。 变 异 过 程 中 的 个 体 有 三 种 情况 ， 如 图 
b) 评 价 。 利 用 评价 函数 评价 中 每 一 个 个 体 。 3 所 示 。 对 于 情况 1， 因 为 不 出 边界 ， 直 接 采 用 式 〈8) 进行 变 
0) 预 选择 机制 。 在 目标 空间 中 使 用 非 支配 解 排 序 和 决策 空 。 异 , 变异 前 个 体 P 变异 成 个 体 ; 对 于 情况 2, 首先 采用 式 (8) 
间 中 使 用 拥挤 距离 选择 方法 ， 选 择 前 50% 的 优良 个 体 并 将 其 放 ”进行 一 次 变异 ， 所 得 变异 个 体 超出 边界 ， 此 时 利用 变异 策略 式 
入 0 中 。 (9) 重新 进行 二 次 变异 ， 新 变异 个 体 没 有 超出 边界 ， 所 以 个 体 
d) 变 异 。 对 8 中 的 每 一 个 个 体 采用 DE/rand12 策略 进行 变 。 ”PP 变异 成 个 体 B 后 再 次 变异 成 个 体 ， 对 于 情况 3， 二 次 变异 
异 。 解 决 变异 边界 问题 的 方法 具体 参见 3.2 节 。 之 后 仍然 超出 边界 ， 将 个 体重 置 为 种 群 内 随机 个 体 ， 因 此 个 体 
6) 交 叉 。 利 用 二 进 制 交叉 算 子 产生 子 种 群 R 。 P 变异 成 个 体 再 次 变异 成 个 体 R， 最 后 设置 为 种 群 中 的 随机 
们 选择 。 在 目标 空间 中 使 用 非 支 配 解 排序 和 拥挤 距离 选择 ”个 体 o 。 
机 制 ， 选 择 w 个 个 体 进入 下 一 代 。 
名 循环 。 如 果 满足 终止 条 件 ， 停 止 并 返回 全 部 最 优 解 ， 否 
则 重新 返回 步骤 o)。 
3.1 个 体 预选 择机 制 
在 变异 之 前 ， 通 过 设计 一 种 个 体 预 选择 机 制 ， 将 初始 种 群 
进行 预 处 理 ， 优 先 挑选 Pareto 最 优 前 沿 和 决策 空间 中 拥挤 距离 
较 大 的 解 ， 这 样 可 以 有 效 增加 所 得 解 的 多 样 性 。 预 选择 机 制 的 . 
过 程 如 图 2 所 示 。 首 先 将 种 群 P 中 每 个 个 体 的 适应 度 值 经 过 非 I 
支配 解 排序 得 出 排名 F1、F2、F3 等 ， 颜 色 越 深 代表 Pareto 前 
沿 排序 越 靠 前 ， 接 着 利用 拥挤 距离 选择 机 制 保留 临界 个 体 并 合 
弃 拥 挤 距 离 小 的 个 体 ， 然 后 得 到 种 群 0 ， 并 利用 种 群 2 进行 变 
异 。 
3.2 ”变异 边界 处 理 方 法 
为 了 提高 种 群 多 样 性 并 保证 足够 的 搜索 范围 ， 本 文采 用 
DE /rand/12 进行 变异 。 其 表达 式 为 
Vo=X, ,+F*(X, —X, )+F*(X, .—-X,,) (9) 种 群 Q 
‘GT Ko #0 Kuo dG HG (情况 3， 出 边界 两 次 ) 
其 中 : Xe、 Xeo、 Xeo、 Xe 和 Xe 是 从 [LN] 中 随机 生成 图 3 变异 边界 的 三 种 情况 和 处 理 方法 
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4 ”实验 及 结果 


4.1 测试 函数 描述 
为 了 测试 MMODE 的 有 效 性 ， 选 择 文献 [12] 中 的 五 个 测试 


函数 进行 测试 ， 分 别 是 F12、F15、F19、F27、F38。 这 五 个 测 
试 函 数 解 的 个 数 如 表 1 所 示 。 这 五 个 测试 函数 的 数学 表达 式 和 


解 的 具体 信息 如 下 : 
表 1 被 测试 的 五 个 非 线性 方程 组 
函数 F12 F15 F19 F27 F38 


问题 1 (F12): 


(10) 


I+1=0 


sin(¥ )-3x1 -1=0 
COs(3X7% )- 


其 中 e[-2.2],1=1.…mn=2。 该 方程 组 有 10 个 解 ， 分 别 是 


(—1.81089, —0.34909) 


(—1.81089, 0.34909) (—1.50222, —0.40908) 


a » 


(-1.50222, 0.40908) (-1.79130,0.30193) (-1.79130,-0.30193) 


9 


(-0.94727,0.78502) (-0.94727,-0.78502) (—0.21306, 1.25685) 


(—0.21306, —1.25685) 。 
问题 2 (F15): 


+ 如一 1.0=0 
台 + 如 一 1.0=0 
+ 区 -1.0=0 
好 二 如 -10=0 
4.731*10 xixs —0.3578x5x — 0.1238x 十 加 
—1.637*10™ x, —0.9338x, -0.3571= 0 
0.2238xixs + 0.7623x,X3 + 0.2638X —N 
一 0.07745 —0.6734x, 一 0.6022 =0 
2 如 十 0.3578x +4.731*10 x, =0 
—0.7623x + 0.2238x, +0.3461= 0 


其 中 : e[-41,i=1.…mn=8 。 该 方程 组 有 16 个 解 ， 分 别 是 
(0.1644,-0.9864,-0.9471 -0.321, -0.9982, -0.0594,0.411,0.9116) ， 


(11) 


(0.1644, —0.9864, —0.9471, -0.321, —0.9982, 0.0594, 0.411, -0.9116) ， 
(0.1644, -0.9864, -0.9471 —0.321,0.9982, —0.0594, 0.411,0.9116) ， 
(0.1644, -0.9864, -0.9471 —0.321,0.9982, 0.0594, 0.411, -0.9116) ， 

(0.1644,—0.9864, 0.7185, —0.6956, —0.998, —0.0638, —0.5278, 0.8494) ， 
(0.1644,—0.9864, 0.7185, —0.6956, —0.998, 0.0638, —0.5278, —0.8494) ， 
(0.1644, —0.9864, 0.7185, —0.6956, 0.998, —0.0638, —0.5278, 0.8494) ， 
(0.1644, —0.9864, 0.7185, —0.6956, 0.998, 0.0638, —0.5278, —0.8494) ， 
(0.6716, 0.741, -0.6516, —0.7586, —0.9625, —0.2711, -0.4376, 0.8992) ， 
(0.6716,0.741, -0.6516, —0.7586, —0.9625, 0.2711, —0.4376, —0.8992) ， 

(0.6716,0.741,—0.6516, —0.7586, 0.9625, —0.2711, -0.4376, 0.8992) ， 

(0.6716,0.741, -0.6516,—0.7586, 0.9625, 0.2711, -0.4376, —0.8992) ， 
(0.6716, 0.741,0.9519, —0.3064, —0.9638, —0.2666, 0.4046, 0.9145) ， 
(0.6716, 0.741,0.9519, —0.3064, —0.9638, 0.2666, 0.4046, -0.9145) ， 
(0.6716, 0.741, 0.9519,-0.3064,0.9638, 0.2666, 0.4046, —0.9145) ， 


(0.6716, 0.741 0.9519, -0.3064, 0.9638, —0.2666, 0.4046, 0.9145) 。 


问题 3 (F19): 
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+O -2=0 12 
D+/4-1=0 (12) 


其 中 ，%e[-22].1=4..mn=2 。 该 方程 组 有 4 个 解 ， 分 别 是 
(~—0.8165,—1.1547), (0.8165,—1.1547) , (-0.8165,1.1547) , (0.8165,1.1547) 。 
问题 4 (F27 ): 
5 一 6 TXT202=0 
—2x6x, + IKON + Xx =0 (13) 
x? +x? —0.265625=0 
其 中 尺 e[-0.6,6],% e[-0.6,0.6],xe[-5,5] 。 该 方程 组 有 12 个 解 ， 分 
别 是 
(0.27986, 0.43279, —0.01419) ， (0.27986, —0.43279, —0.01419) ， 
(-0.27986, 0.43279, —0.01419) , (-0.27986,-0.43279,-0.01419) ， 
(0.46698, 0.21807, 0.00000) , (—0.46698, 0.21807, 0.00000) ， 
(0.46698, —0.21807, 0.00000) , (—0.46698, —0.21807, 0.00000) ， 
(0.00000, 0.51539, 0.00000) ， (0.00000, =—0.51539, 0.00000)， 
(0.51539, 0.00000, -0.01245) , (-0.51539, 0.00000, -0.01245) 。 
问题 5 (F38): 
oy +4m -6.0=0 
-0.6787 =0 (4) 
其 中 :ae[-23,es[QL]] 。 该 方程 组 有 4 个 解 ， 分 别 是 


(—1.56353, 0.27763) (—0.78971,1.08830) (1.56353,0.27763) 


% 


(0.78971,1.08830) 。 
4.2 评价 指标 

对 于 有 限 解 的 非 线 性 方程 组 ， 需 要 在 算法 运行 结束 后 确定 
己 搜 索 到 解 的 个 数 。 以 其 中 一 个 解 为 例 ， 如 果 算 法 获得 的 解 与 
已 知 最 优 解 的 欧 氏 距离 小 于 预先 设 定 的 精确 度 * ， 就 认为 该 最 
优 解 已 经 找到 。 这 里 将 精确 度 设 为 


时 0.01 nz<5 
“2101 其 他 (15) 


其 中 : 
的 个 数 之 间 的 比值 称 为 峰 比率 ， 


中 找到 最 优 解 的 平均 概率 。 如 果 全 部 最 优 解 都 能 在 一 次 运行 中 
找到 ， 那 么 就 称 为 一 次 成 功 的 运行 。 成 功 运行 的 次 数 与 全 部 运 


站 


了 总 次 数 之 间 的 比值 称 为 成 功率 ， 用 SR 表示 ， SR 代表 单 次 运 
行 成 功 找到 全 部 最 优 解 的 平均 概率 。 本 文采 用 PR 和 SR 两 个 性 
能 指标 来 评价 当前 算法 的 性 能 ， 它 们 的 值 越 大 ， 就 说 明 算 法 的 
性 能 越 好 。 
4.3 比较 算法 

与 本 文 所 提出 的 MMODE 比较 的 算法 如 下 : 

a)A-WeB。 该 算法 由 Gong 等 人 “于 2017 年 提出 ， 主 要 是 
多 权 值 双 目标 转换 技术 求解 非 线性 方程 组 的 最 优 解 。 

b)DN-NSGAI。 该 算法 由 Yue 等 人 "” 于 2016 年 提出 , 主要 
是 用 决策 空间 非 支配 解 排序 算法 求解 多 模 态 多 目标 优化 问题 。 

cj)NSGAII。 该 算法 由 Deb 等 人 "于 2002 年 提出 ， 用 于 求 

解 非 支 配 解 排序 的 问题 。 
d)Omni-Opt。 该 算法 1 


Deb 等 人 "于 2005 年 提出 , 用 于 改 
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善 搜索 解 在 决策 空间 的 分 布 。 
4.4 参数 设置 
在 实验 中 参数 设置 如 下 : 种 群 大 小 NP=100, 最 大 评价 次 数 
max FES =50000, 缩放 因子 F =0.5, 交叉 率 Cr =0.9。 为 保证 比较 
的 公平 性 ， 所 有 算法 的 种 群 大 小 、 最 大 评价 次 数 等 参数 均 保持 
相等 。 每 种 算法 均 独 立 运 行 51 次 ， 结 果 取 平均 值 。 
4.5 实验 结果 和 分 析 

为 避免 偶然 性 ， 首 先 将 第 26 次 所 提 MMODE 算法 搜索 到 
的 非 线性 方程 组 的 解 的 情况 输出 ， 其 中 问题 1 一 5 的 所 有 解 在 
决策 空间 和 目标 空间 的 分 布 如 图 4 所 示 。 从 图 4(a) 可 以 看 出 ， 
对 于 有 2 个 决策 变量 的 问题 1， 在 决策 空间 和 目标 空间 上 ，10 
个 解 都 能 够 找到 , 而 且 搜索 到 的 解 与 真实 最 优 解 之 间 比 较 接 近 ， 
可 以 看 出 MMODE 对 于 问题 1 解决 的 比较 好 ， 它 的 10 个 最 优 
解 全 部 都 能 搜索 到 ; 从 (b) 可 以 看 出 ， 对 于 有 8 个 决策 变量 的 
问题 2， 由 于 维度 的 限制 ， 这 里 给 出 了 前 两 维 的 图 形 ， 在 决策 
空间 上 , 搜索 到 的 解 都 分 布 在 16 个 最 优 解 的 周围 , 而 在 目标 空 
间 上 ， 搜 索 到 的 解 与 真实 最 优 解 之 间 还 有 一 定 的 差距 ， 可 以 看 
出 MMODE 对 于 问题 2 解决 的 并 不 好 ， 它 的 16 个 最 优 解 还 不 
能 很 好 的 搜索 到 , 说 明 MMODE 对 于 决策 变量 较 多 的 问题 2 的 
求解 上 仍然 还 有 竺 改进， 从 (c) 可 以 看 出 ， 对 于 有 2 个 决策 变 
量 的 问题 3, 在 决策 空间 和 目标 空间 上 , 4 个 解 都 能 够 找到 ， 而 
且 搜 索 到 的 解 与 真实 最 优 解 之 间 比 较 接近 ， 不 过 仍 存在 一 些 与 
真实 最 优 解 相 距 较 远 的 搜索 解 ， 可 以 看 出 MMODE 对 于 问题 3 
解决 的 比较 好 ， 它 的 4 个 最 优 解 全 部 都 能 搜索 到 ， 但 是 搜索 到 
的 全 部 解 还 不 能 完全 地 趋 近 真实 最 优 解 ， 从 Cd) 可 以 看 出 ， 对 
于 有 3 个 决策 变量 的 问题 4， 在 决策 空间 和 目标 空间 上 ，12 个 
真实 最 优 解 只 有 8 个 能 够 被 找到 ， 另 外 4 个 最 优 解 与 搜索 到 的 
解 还 有 一 些 差距 , 可 以 看 出 MMODE 对 于 问题 4 只 能 搜索 到 一 
部 分 真实 最 优 解 ， 从 (e) 可 以 看 出 ， 对 于 有 2 个 决策 变量 的 问 
题 5, 在 决策 空间 和 目标 空间 上 , 4 个 真实 最 优 解 都 能 被 很 好 地 
找到 ， 可 以 看 出 MMODE 对 于 问题 $ 的 求解 表现 的 比较 好 。 

接 下 来 , 将 MMODE 算法 与 A-WeB 、DN-NSGAII NSGAII 
和 Omni-Opt 算法 在 峰值 率 上 进行 比较 , 详细 结果 如 表 2 所 示 。 
从 表 2 可 以 看 出 ， 对 于 问题 1，MMODE 在 峰 比 率 PR 上 数值 最 
大 ，A-WeB 排 第 二 ， A-WeB 和 MMODE 相差 不 大 ,NSGAII、 
Omni-Opt 和 DN-NSGAII 分 别 排 在 第 三 到 第 五 ， 说 明 MMODE 
平均 搜索 到 解 的 个 数 比 其 他 算法 更 多 ; 对 于 问题 2, MMODE 在 
峰 比 率 PR 上 数值 最 大 , A-WeB 排 在 第 二 , NSGAIIDN-NSGAII 
和 Omni-Opt 分 别 排 在 第 三 到 第 五 ， 与 前 两 名 相差 较 多 ， 说 明 
MMODE 平均 搜索 到 的 解 的 个 数 比 其 他 算法 更 多 ， 并 且 效 果 比 
其 他 算法 表现 更 好 ; 对 于 问题 3，MMODE 在 峰 比 率 PR 上 数值 
为 1，A-WeB、NSGAIT、DN-NSGAII 和 Omni-Opt 分 别 排 在 第 
二 到 第 五 ， 但 与 MMODE 相差 不 大 ， 说 明 MMODE 每 次 都 能 
够 搜索 到 全 部 最 优 解 ， 其 他 四 种 算法 也 可 以 搜索 到 绝 大 部 分 的 
最 优 解 ; 对 于 问题 4, MMODE 在 峰 比 率 PR 上 数值 最 大 , Omni- 
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说 明 MMODE 在 五 个 算法 比较 上 表现 最 好 , 但 是 全 部 算法 平均 
搜索 到 的 最 优 解 并 不 多 ; 对 于 问题 5, MMODE 和 A-WeB 在 峰 
比率 PR 上 数值 均 为 1, NSGAII 排 在 中 间 , DN-NSGAI 和 Omni- 
Opt 的 PR 值 一 样 大 ， 排 在 最 后 ， 说 明 MMODE 和 A-WeB 每 一 
次 都 能 够 搜索 到 全 部 最 优 解 ， 其 他 三 个 算法 可 以 搜索 到 大 部 分 


的 最 优 解 。 


问题 1 PS 


问题 1 PF 


1.5r 2.5 
1 | ， 真实 最 优 解 中 过 *| 真实 最 优 解 
[搜索 的 解 ， | 搜索 到 gj 和解 
四 二 1.5 
No 
四 1 
-0.5 
引 0.5 
引 生 15 区 | -05 0 05 0.5 0 0.5 1 
测 f1 
(a) 问 题 1 的 10 个 最 优 解 在 决策 空间 和 目标 空间 的 分 布 
问题 2 PS 问题 2 PF 
^ 搜索 到 的 解 1.5 
0 a 二 
1 二 真实 最 优 解 
中 搜索 到 的 解 
0 02 04 06 08 1 -05 05 


0 
f1 


(b) 问 题 2 的 16 个 最 优 解 如 


问题 3 PS 
2 


FE 决策 气 间 和 


some 


SN 0 丰富 
x 0 。 搜索 到 的 


人 0 0 05 
对 


RN1 


问题 3 PF 


1 


标 空 间 的 分 布 


mom ™ 


”真实 最 优 解 


^ 搜索 到 的 解 


1 4 


-0.5 


0 
fl 


0.5 


1 


(0) 问 题 3 的 4 个 最 优 解 有 有 


FE 决 策 空 间 和 


标 空间 的 分 布 


a4 PS 间 题 4 PF 
15 
5 Re 真实 最 优 解 
+| 真实 最 优 解 搜索 网 角 
0 ws. “所 寺 玖 争 | & 1 本 
5 ny 
05 
0 2 4 0 和 5 | 
D4 -02 0 02 04 
5 0 1 


(gd) 问 题 4 的 12 个 最 优 解 如 


标 空 间 的 分 布 


ja PS 问题 5 PF 
| 2 
1 二 真实 最 优 解 | 真实 最 优 解 
^ 搜索 到 的 解 4 | 搜索 到 的 解 
Nos5 a1 g 
O01 
x f 
(e) 问 题 5 的 4 个 最 优 解 在 决策 空间 和 目标 空间 的 分 布 


表 2 MMODE 和 


图 4 问题 1~5 的 所 有 人 解 在 决策 空间 不 


瑟 


标 空间 的 分 布 


个 算法 在 问题 1-5 上 的 峰 比 率 PR 


算法 问题 1 问题 2 问题 3 问题 4 问题 5 
MMODE 0.8980 0.7350 1.0000 0.3967 1.0000 
A-WeB 0.8880 0.6688 0.9950 “0.0933 1.0000 
DN-NSGAII 0.4080 0.0850 0.9650 0.2800 0.7100 
NSGAII 0.6400 0.0888 0.9800 0.2783 0.7250 
Omni-Opt 0.5020 0.0725 0.9550 0.2900 0.7100 


Opt、DN-NSGAII、 NSGAII 和 A-WeB 和 分 别 排 在 第 二 到 第 五 ， 
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为 了 避免 单独 比较 峰值 率 造成 的 偏差 ， 现 在 将 五 种 算法 在 
成 功率 上 进行 比较 ， 详 细 结 果 如 表 3 所 示 。 从 表 3 可 以 看 出 ， 
对 于 问题 1，MMODE 和 A-WeB 在 成 功率 SR 上 数值 一 样 ， 排 
在 第 一 ，DN-NSGAII、NSGAII 和 Omni-Opt 数值 为 零 ， 说 明 
MMODE 和 A-WeB 成 功 找到 全 部 最 优 解 的 次 数 最 多 ， 但 是 成 
功率 并 不 高 ， 对 于 问题 2，MMODE 排 在 第 一 ， 其 他 四 种 算法 
成 功率 均 为 零 ， 说 明 MMODE 可 以 找到 全 部 最 优 解 , 但 是 成 功 
率 非常 低 ; 对 于 问题 3, MMODE 成 功率 是 1, A-WeB、NSGAII、 
DN-NSGAI 和 Omni-Opt 分 别 排 在 第 二 到 第 五 ， 说 明 MMODE 
每 次 运行 都 可 以 成 功 找到 全 部 最 优 解 ， 对 于 问题 4，A-WeB 成 
功率 排 在 第 一 , MMODE、DN-NSGAII、NSGAII 和 Omni-Opt 
成 功率 均 为 零 , 说 明 MMODE 每 次 运行 都 不 能 成 功 找到 全 部 最 
优 解 ， 对 于 问题 5, MMODE 和 A-WeB 排 在 第 一 ， 并 且 成 功率 
都 是 1，NSGAII、DN-NSGAI 和 Omni-Opt 依次 排 在 后 三 位 ， 
与 前 两 名 相差 较 大 ， 说 明 MMODE 和 A-WeB 每 次 运行 都 能 成 
功 找 到 全 部 最 优 解 。 


表 3 MMODE 和 四 个 算法 在 问题 1-5 上 的 成 功率 SR 


算法 问题 1 问题 2 问题 3 问题 4 问题 5 
MMODE 0.2800 0.0400 1.0000 0.0000 1.0000 
A-WeB 0.2800 0.0000 0.9800 0.0933 1.0000 
DN-NSGAI 0.0000 0.0000 0.8600 0.0000 0.0600 
NSGAII 0.0000 0.0000 0.9200 0.0000 0.0800 
Omni-Opt 0.0000 0.0000 0.8200 0.0000 0.0400 


综合 表 2 和 3 可 以 看 出 ， 对 于 大 部 分 测试 问题 , MMODE 
在 峰 比 率 和 成 功率 上 的 结果 比 其 他 四 种 算法 表现 更 好 ， 说 明 
MMODE 在 求解 非 线性 方程 组 问题 上 是 有 效 的 。 

另外 ， 为 了 更 好 地 观察 实验 数据 的 分 布 ， 对 其 中 的 四 种 算 
法 的 峰值 率 PR 所 对 应 的 盒 形 图 表示 如 图 5 所 示 。 其 中 横 坐 标 
1、2、3、4 分 别 代表 MMODE、DN-NSGAII、NSGAII、Omni- 
Opt 四 种 算法 ， 纵 坐标 代表 峰值 率 PR 值 。 从 图 5 可 以 看 出 ， 对 
于 问题 1 一 4，MMODE 和 其 他 三 种 算法 的 波动 范围 大 致 相同 ， 
但 峰值 率 的 均值 比 其 他 三 种 算法 更 大 ， 对 于 问题 3， 四 种 算法 
的 峰值 率 波动 性 比较 小 ， 峰 值 率 都 接近 于 1， 但 是 MMODE 比 
其 他 三 种 算法 更 加 稳定 ， 每 次 都 可 以 找到 全 部 最 优 解 ， 对 于 问 
题 5， 四 种 算法 的 峰值 率 波 动 性 也 比较 小 ， 但 MMODE 的 峰值 
率 可 以 一 直 保持 为 1， 这 比 其 他 三 种 算法 长 期 保持 在 0.7 附近 
效果 更 好 。 因 此 ， 从 盒 形 图 上 看 ，MMODE 在 大 部 分 问题 上 比 
其 他 算法 表现 更 好 。 
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MMODE 的 显著 性 优势 从 小 到 大 依次 是 A-WeB、NSGAI、 
Omni-Opt 和 DN-NSGAIL 说 明 MMODE 算法 显著 性 地 优 于 其 
他 四 种 算法 。 
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图 5 四 种 算法 的 盒 形 图 


表 4 MMODE 和 四 个 算法 在 问题 1~5 上 的 Wilcoxon 秩 和 检验 结果 


MMODE vs P 值 
A-WeB 3.7429e-13 
DN-NSGAII 6.0879e-25 
NSGAII 4.4421e-19 
Omni-Opt 4.4891e-24 
5 应 用 
非 线 性 方程 组 在 工程 中 也 有 广泛 的 应 用 ， 比 如 可 以 用 来 建 


立 两 个 连续 的 非 绝缘 搅拌 等 式 反 应 器 的 数学 模型 ， 从 而 得 出 系 
统 稳定 后 两 个 反应 器 的 温度 “。 这 些 反应 往往 处 在 一 个 稳定 可 


为 了 进一步 比较 MMODE 算法 对 其 他 四 种 算法 的 优势 , 这 
里 使 用 Wilcoxon 秩 和 检验 进行 分 析 。 因为 考虑 到 搜索 到 的 非 线 
性 方程 组 解 的 平均 数量 是 当前 判断 算法 好 坏 的 重要 标志 ， 所 以 
这 里 检验 的 对 象 是 MMODE 和 其 他 四 种 算法 的 全 部 51 次 峰值 
率 。 通 过 进行 Wilcoxon 秩 和 检验 ， 比 较 结 果 如 表 4 所 示 ， 其 中 
全 部 检验 结果 都 是 在 显著 性 水 平 为 0.05 下 进行 的 。 从 表 4 可 以 
看 出 ,与 MMODE 比较 的 四 种 算法 得 到 的 P 值 都 远 远 小 于 0.05， 


循环 的 状态 ， 并 且 是 一 个 一 阶 不 可 道 的 放 热 反应 。 忽 略 个 别 因 
素 的 影响 ， 该 反应 器 的 数学 模型 构建 如 下 : 
Dp 10x 
fe 


(16) 


D 10x, 
f=n +p) +10- 用 | 名 -ps 二 人 Bs ep 


其 中 : 石和 浆 分 别 代表 两 个 反应 器 中 的 无 量 纲 温度 ; /表示 无 
量 纲 激 活 能 量 ，y=1000 ; D 表示 达 姆 科勒 数 ，D=22; hh 和 忆 
分 别 表 示 无 量 纲 传 热 系数 ， B=P=2;R 表示 循环 率 ， R=0.96 。 
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现 将 MMODE 和 四 种 算法 分 别 用 来 求解 这 一 工程 实例 , 讨 
价 指标 峰 比 率 和 成 功率 的 结果 如 表 5 所 示 。 
表 5 MMODE 和 四 个 算法 在 工程 实例 上 的 结果 


评价 指标 MMODE A-WeB DN-NSGAIT NSGAI Omni-Opt 
峰 比 率 0.8571 0.8371 0.8367 0.8385 0.8543 
成 功率 0.1569 0.1176 0.0784 0.0784 0.0980 


从 表 5 可 以 看 出 ， 在 峰 比 率 上 ， 五 种 算法 的 结果 总 体 相 差 
较 小 ，MMODE 峰 比 率 最 大 ; 在 成 功率 上 ， 五 种 算法 的 成 功率 
都 不 高 ， 但 相 比 较 而 言 ，MMODE 有 最 大 的 成 功率 。 综 合 比 较 
条 比率 和 成 功率 的 结果 可 以 得 出 ，MMODE 在 非 绝 缘 搅 拌和 釜 式 
反应 器 上 的 应 用 也 比 其 他 四 种 算法 具有 更 好 的 有 效 性 。 
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本 文 主要 介绍 了 非 线性 方程 组 的 求解 ， 首 先 总 结 了 将 非 线 
生 方 程 组 转换 为 优化 问题 的 主要 方法 ， 然 后 根据 当前 算法 存在 
的 问题 ， 提 出 了 一 种 多 模 态 多 目标 差分 进化 算法 。 该 算法 的 主 
要 改进 包括 种 群 预选 择机 制 和 边界 选取 方法 ， 这 些 改进 能 够 增 
加 决策 空间 解 的 多 样 性 。 为 了 验证 多 模 态 多 目标 差分 进化 算法 
的 有 效 性 , 对 5 个 测试 函数 进行 了 51 次 独立 实验 。 最 后 经 过 与 
四 种 算法 得 出 的 结果 进行 比较 ， 可 以 说 明 多 模 态 多 目标 差分 进 
化 算法 在 峰 比 率 和 成 功率 两 个 指标 上 都 有 一 定 的 优势 。 另 外 ， 
通过 求解 一 个 具体 的 工程 实例 ， 也 进一步 验证 了 所 提 算 法 的 有 
效 性 。 因 此 ， 多 模 态 多 目标 差分 进化 算法 是 一 种 有 效 求解 非 线 
性 方程 组 的 方法 。 不 过 ， 对 于 非 线性 方程 组 问题 的 求解 ， 尤 其 
是 在 成 功率 上 ， 当 前 算法 仍然 还 有 很 大 的 提升 空间 。 下 一 步 工 


Se 


作 将 继续 优化 算法 ， 使 之 能 够 更 好 地 应 用 到 更 多 复杂 的 工程 实 
例 中 。 
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